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METHOD AND SYSTEM FOR REAL 
USING SATELLITE TRANSMITTED THREE-CARRIER RADIO 
SIGNALS AND IONOSPHERIC CORRECTIONS 

[0001] The invention concerns a method and a system for real-time navigation 
using satelHte-transmitted three-carrier radio signals and ionospheric corrections, more 
precisely corrections obtained by means of a continuously updated real-time ionospheric 
model, the model being based on data from a satellite navigation system, for example 
implemented like a three-dimensional voxel model. 

[0002] It is particularly though not exclusively, applicable to the field of high- 
precision instantaneous navigation, typically having a precision within one decimeter, sls 
will be shown below, at distances on the order of hundreds of kilometers or more. 

[0003] By way of illustration, we will hereinafter focus on the preferred 
application of the invention, without in any way limiting its scope. 

[0004] One of the current commonly used techniques for obtaining a precise 
positioning of an "object," whether stationary or moving, and in the latter case, data 
relative to its position, its movement, the direction of this movement and/or its speed, is to 
use radio signals transmitted by artificial satellites orbiting around the earth. The term 
"object" should be understood in its most general sense, notably a land, sea or air vehicle. 
For simplicity's sake, we will hereinafter refer to this "object" as a "rover." 

[0005] There are various known techniques for obtaining the aforementioned 
positioning. They are based in particular on the knowledge of the instantaneous position 
of several satellites in space (or a constellation of satellites, as these satellites may or may 
not be geostationary) and of the propagation speed of the radio waves. High-precision 
clocks are installed aboard the satellites, and the transmitted signals include time-stamped 
information, which makes it possible to know precisely the moment of both transmission 
and reception. It is thus possible to determine the theoretical distance separating a 
satellite in view of the rover from the latter at a given moment, knowing the propagation 
speed of the waves and the time they take to reach the rover. If a sufficient number of 
satellites is observable, it is thus possible to determine the coordinates of the rover 
relative to a reference frame, in two dimensions (longitude and latitude on the earth) or 
even three dimensions (longitude, latitude and altitude / vertical). 

[0006] However, as will be shown, because of the accumulation of errors in the 
measurements due to various causes, the distances calculated are only approximate, and 
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the determination of the position of the rover suffers from a more or less substantial 
degree of imprecision, depending on the technologies used. 

[0007] One of the best known satellite navigation systems is the system known as 
"GPS" for "Global Positioning System." 

[0008] Customarily, "GPS" (or similar system) satellites transmit in two 
frequency bands, generally designated Li (a carrier frequency equal to 1.575 GHz) and L2 
(a carrier frequency equal to 1.227 GHz), hence the adjective "dual frequency" that is 
applied to them. 

[0009] The use of these two frequencies, in accordance with certain methods well 
known to one skilled in the art, makes it possible to improve the precision of the 
determination of a rover's position relative to a reference frame, but it requires "GPS" 
receivers that are more complex and more costly. 

[0010] A position determination can be made using two main methods: real-time, 
or after-the-fact, by performing what is known as a "post-processing." The first case is 
commonly referred to as "single-epoch" resolution (a term that will be used below) or 
instantaneous resolution, the calculations being performed during a single observation 
"epoch." The second method ("post-processing") makes it possible to improve precision. 
However, while the latter method does not present any major drawbacks for slow-moving 
rovers (ships, for example), it is not appropriate for rovers that move very fast (aircraft, 
for example). 

[001 1] Precision can be further improved by combining the signals transmitted by 
the satellites with signals originating from fixed ground reference stations whose 
positions are well known. However, if the rover travels great distances, it is necessary for 
this network of stations to be relatively dense, especially in cases where high precision is 
desirable in the determination of the rover's position, which correspondingly increases the 
cost of the global system. 

[0012] Moreover, among the numerous causes of errors, differential ionospheric 
refraction, when considering distances equal to tens of kilometers or more, is one of the 
main problems affecting capacities for instantaneous resolution of carrier phase 
ambiguity, and consequently the capability to provide a navigation wherein the precision 
is on the order of one centimeter with dual-frequency global navigation satellite systems 
such as the aforementioned "GPS" system. This characteristic will remain true with 
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respect to future three-frequency systems like the "GALILEO" system and the 
"Modernized GPS" system. 

[0013] In essence, the three-carrier systems currently being planned offer the 
potential advantages of a high success rate and high integrity in instantaneous ambiguity 
resolution, with a minimal number of geodetic calculations. This is particularly due to the 
fact that a higher quantity of different data (i.e., associated with the aforementioned three 
frequencies) is made available, which correspondingly improves the chances of obtaining 
an instantaneous ("single epoch") ambiguity resolution. 

[0014] But here again, this resolution can be seriously affected by ionospheric 
refraction, as explained below. 

[0015] In order to achieve high precision in the instantaneous determination of the 
position of a rover, particularly a rover that travels great distances, there is still a need to 
implement techniques that make it possible, in particular, to reduce the harmful influence 
of ionospheric refraction. 

[0016] In the prior art, various methods for meeting this need have been proposed. 

[0017] For example, there is the method known as "TCAR," for "three-carrier 
ambiguity resolution." This method is described in the article entitled "ANALYSIS OF 
THREE-CARRIER AMBIGUITY RESOLUTION (TCAR) TECHNIQUE FOR 
PRECISE RELATIVE POSITIONING IN GNSS-2," by U. VOLLATH et al., published 
in "Proceedings of the ION GPS" 1998, IX-O-13, pages 1-6. 

[0018] There is also the method known as "CIR," for "cascade integer 
resolution." This method is described in the article by Jaewo JUNG et al., entitled 
"Optimization of Cascade Integer Resolution with Three Civil GPS Frequencies," 
published in "Proceedings of the ION GPS 2000." 

[0019] These two techniques share a similar basic approach: the double difference 
ambiguities of integers are successively resolved by calculating wave frequency beats. 
This calculation is performed from the longest to the shortest beat wavelength, including 
combinations of so-called "wide" lane and "extra wide" lane carrier phases (with 
wavelengths of 7.480 m and 0.862 m, respectively), and a first carrier at the "Li 
frequency" (with a wavelength of 0.190 m). 

[0020] The "TCAR" method in particular constitutes a simple approach that tries 
to resolve the frill set of ambiguities instantaneously (in "single epoch" mode). But the 
performance of "TCAR" is strongly affected by the ionospheric refraction decorrelation 
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that occurs with distance. In fact, as explained below, ionospheric delay is a problem 
when (as in the case of two-frequency systems) the value of its double differential is more 
than 0.26 TECU (which corresponds to a 4-cm delay for Li). 

[0021] A "TECU" is a unit used to describe certain electrical characteristics of the 
ionosphere. In essence, the ionosphere can be described using a map that represents a 
count of the total number of electrons, or "TEC" (for "Total Electron Content"). The map 
represents the integration of the number of electrons in a vertical direction as a function 
of latitude and longitude. A unit of TEC is referred to as a "TECU" (for "TEC Unit"), 
with one TECU = lO'^ electrons contained in a cylinder aligned on the line of observation 
of an observed satellite, the cross-section of which is 1 m^. The charged particles in the 
ionosphere are generated by the sun, whose radiation intensity varies naturally as a 
function of the time in question. Because the earth rotates on its axis underneath the 
ionospheric layer, the "TEC" map is normally considered to represent a reference frame 
that is fixed relative to the sun, but that changes as a function of time. 

[0022] The above-mentioned threshold is easily exceeded, as may be seen by 
consulting ionospheric ("TEC") maps of vertical delays calculated from "GPS" data. 
Such maps are issued, for example, by the "Jet Propulsion Laboratory," the University of 
Bern, etc., and published on the Intemet by the "University Corporation for Atmospheric 
Research" and other similar organizations. 

[0023] Consequently, in order to further improve the "TCAR" method, an 
integrated approach known as "ITCAR" (for "Integrated TCAR") was developed. This 
technique is described, for example, in the above-mentioned article by VOLLATH et al. 

[0024] This technique uses search algorithms and a navigation filter wherein the 
ambiguities are part of the output signals and the residual ionospheric errors are roughly 
estimated. For a more detailed description of the techniques used, it would be worthwhile 
to refer to this article. 

[0025] However, although it provides a significant improvement, the "ITCAR" 
technique is nonetheless still affected by the lack of knowledge of the double difference 
of the ionospheric refraction, thus limiting the success rate of the ambiguity resolution for 
distances greater than several tens of kilometers, as described in the article by VOLLATH 
et al., entitled "Ambiguity Resolution Using Three Carriers - Performance Analyzing 
Using 'Real' Data, published in "GNSS Symposium," Seville, May 2001. 
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[0026] It has also been proposed, again with a view to improving the precision of 
the determination of a rover's position relative to a reference frame, to combine a real- 
time ionospheric model of the ionosphere, obtained from "dual-frequency" data generated 
by a network of fixed stations, with data from a geodetic program, and to use such data to 
perform ionospheric corrections. This method has been used with some success in 
resolving ambiguities in real-time in two-frequency systems of the "GPS" type. 

[0027] One method of this type, called "WARTK" (for "Wide Area Real-Time 
Kinematics"), is described for example in the article by Hemandez-Pajares et al. entitled 
"Tomographic Modeling of GNSS Ionospheric Corrections: Assessment and Real-Time 
Applications," published in "ION GPS" 19-22 September 2000, pages 616-625. This 
method does make it possible to attenuate the harmful effects of the disturbances 
experienced by the radio waves propagating in the ionosphere, and consequently, to 
substantially improve the success rate of the phase ambiguity resolution and the 
determination of a rover's position relative to a reference frame, but it requires a large 
number of calculations to produce the aforementioned model in real time and to 
determine, also in real time, the ionospheric corrections to be applied to the distance 
measurements. 

[0028] The invention seeks to overcome the drawbacks of the methods and 
systems of the prior art, some of which have been mentioned. 

[0029] The object of the invention is to provide a method for high-precision real- 
time navigation (a determination of the position of a rover typically better than one 
decimeter) using three-carrier radio signals transmitted by satellites and ionospheric 
corrections obtained from data in an ionospheric model describing a region passed 
through by the radio radiation, combined with geodetic data. 

[0030] To this end, according to a first important characteristic of the invention, a 
simple method for the instantaneous ("single epoch") resolution of the ambiguities of 
satellite-transmitted three-carrier phases, of the aforementioned "TCAR" type, is used. 
Precise real-time ionospheric corrections, based on a descriptive ionospheric model, are 
generated by fixed reference stations. These ionospheric corrections are transmitted to 
users that incorporate them into "TCAR." 

[0031] The invention has three important advantages. 

[0032] The capabilities specific to the invention make it possible to obtain better 
performance than the "ITCAR" method over long distances, at a position very far from 
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the nearest reference station. The method according to the invention typically allows the 
instantaneous resolution of more than 90% of ambiguities at more than 100 kilometers 
from this station, and practically 100% at a distance on the order of 60 km, even in an 
operating scenario that is not very favorable: low ionospheric values and maximum 
sunlight conditions. Under comparable conditions, the methods of the prior art using 
three frequencies achieve a typical success rate that is only on the order of 60%. 

[0033] It is easy to see that this latter characteristic is very advantageous, since it 
does not require the installation of a very dense network of fixed ground reference 
stations (stations on the order of a hundred kilometers apart). The result is a very 
substantial overall savings, for example in the case of Europe: in covering that continent, 
it is possible to reduce the number of reference stations by 99% compared to a network 
wherein the reference stations are only separated by a distance on the order of ten 
kilometers, while retaining the same precision in the determination of a rover's position 
as that obtained by the methods of the prior art. This is in fact the installation density 
required to obtain a navigation precision better than one decimeter, as sought by method 
of the invention in its preferred application. Moreover, the network of fixed ground 
reference stations required by the method of the invention is very similar to the network 
currently being installed for the European satellite navigation system known by the 
abbreviation "EGNOS GNSS" (for "European Geostationary Navigation Overlay System 
- Global Navigation Satellite Systems"). This system is based on a constellation of 
'TMMERSAT III" satellites, which suggests that there is substantial compatibility with 
the installations that currently exist or are in the process of being installed. 

[0034] Nevertheless, the method according to the invention is still quite simple, 
like the "TCAR" method. In particular, it requires only a small number of calculations 
for the user (the rover's receiver) compared to the "ITCAR" and "WARTK" methods, 
which were developed for two-frequency positioning systems, as explained above. 

[0035] Furthermore, the descriptive real-time ionospheric model of the region 
passed through by the radio radiation transmitted by the observed satellites uses only 
carrier phase data, this data being combined with geodetic estimates calculated at a single 
center, which can be one of the stations, in the reference station network, called a 
"master" station. 

[0036] Hence, the primary subject of the invention is a method for real-time 
navigation using three-carrier radio signals of first, second and third different frequencies 
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that increase in value from said first to said third fi-equency, to determine the position of a 
user, called a rover, said radio signals being transmitted by a given number of transmitters 
installed aboard satellites orbiting around the earth and in view of said rover, said signals 
being received by a receiver associated with said rover and a receiver associated with at 
least one ground station among a plurality of fixed ground stations called reference 
stations, said radio signals passing through a so-called ionospheric layer of the 
atmosphere surrounding said earth and experiencing disturbances that generate phase 
ambiguities in said carriers, characterized in that it comprises at least the following steps: 

[0037] - a first step consisting in the determination, in said rover, of the so-called 
"extra-wide lane" ambiguity of the phase differences between said third and second 
carriers, from a combination of pseudoranges using a single code value; 

[0038] - a second step consisting in the estimation, in said rover, of the so-called 
"wide lane" ambiguity of the phase differences between said first and second carriers, 
from said "extra-wide lane" ambiguity determined during said first step; 

[0039] - a third step consisting, in said rover, in the resolution of the ambiguity of 
one of said firequencies, from said "wide lane" ambiguity estimated during said second 
step; and 

[0040] an additional step for applying real-time ionospheric corrections during 
said third step, said ionospheric corrections being based on a continuously updated real- 
time ionospheric model of said layer. 

[0041] Another object of the invention is a satellite navigation system that 
implements this method. 

[0042] The invention will now be described in greater detail in reference to the 
attached drawings, among which: 

[0043] - Fig. 1 A schematically illustrates a partial cross-section of the earth, and 
of the ionospheric layer that surrounds it, cut into voxels, as well as radio beams 
generated by "GPS" transmitters on three communication satellites, and three ground 
stations receiving these signals; 

[0044] - Fig. IB schematically illustrates the architecture of a complete 
navigation system for implementing the method of the invention, seen fi-om a "GPS" 
receiver of a rover; 

[0045] - Fig. 2 schematically represents the baseline and the positions of a ground 
reference station and a rover, relative to the same coordinate axes. 
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[0046] - Fig. 3 is an enlargement (zoom) of Fig. 2 showing the trajectory of the 
rover relative to these same coordinate axes; 

[0047] - Fig. 4. is a graph illustrating examples of delays estimated from a real- 
time ionospheric model according to the invention compared to real ionospheric delays; 

[0048] - Fig. 5 is a graph illustrating examples of double differences of "Slant 
Total Electron Content" ("STEC") estimated in real time with a descriptive model of the 
ionosphere, compared to real values; 

[0049] - Fig. 6 is a graph illustrating examples of errors in the determination of 
the real-time ionospheric model of the double difference of "STEC" compared to two 
ionospheric thresholds; 

[0050] - Figs. 7A through 7D are graphs illustrating examples of estimated errors 
(in meters) in different navigation components; 

[0051] - Figs. 8 A through 8D are graphs illustrating examples of real errors (in 
meters) in the same navigation components; 

[0052] - Fig. 9 is a graph illustrating an exemplary success rate for the 
instantaneous determination of an ionospheric double difference that is high enough to 
resolve all of the ambiguities (with a "TECU" better than 0.26) for corrections with 
latency times within a range of 1 to 30 seconds; and 

[0053] - Fig. 10 is a graph illustrating an exemplary success rate in the 
instantaneous determination of an ionospheric double difference for corrections with 
latency times of up to 900 seconds. 

[0054] Hereinafter, without in any way limiting the scope of the invention, we 
will focus on the context of its preferred application unless otherwise indicated, i.e., on 
the case of a high-precision satellite navigation method using three-carrier radio waves of 
different frequencies transmitted by satellites, and real-time correictions derived from a 
descriptive real-time ionospheric model of a region sounded by radio frequency radiation 
from satellites. 

[0055] We will now describe the main steps of the method according to the 
invention in detail, using typical numeric values to illustrate them. It comprises three 
basic steps, which are similar to those of the aforementioned "TCAR" method and are 
explained below: 

Step 1 : 
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[0056] In order to resolve the so-called "extra wide lane" ambiguity (with a 
typical wave-length of 7.5 meters in a data set of a given exemplary experiment), a 
combination of codes known as "pseudoranges" is added. 

[0057] Among the known causes of errors is the so-called "multipath" 
phenomenon. Although pseudorange multipaths can reduce the success rate, this error is 
generally of low amplitude compared to the long wavelength of the "extra-wide lane," 
and it is usually possible to overcome this problem. 

[0058] In fact, the first step typical of the "TCAR" approach (see for example the 
above-mentioned article by VOLLATH et al., from 1998) is to estimate the double 
difference of the ambiguity (between satellite-receiver pairs), hereinafter called VAN^^ , 

from the double difference of "extra-wide" lane carrier phases, hereinafter called VAL^^^, 
(Lew being the wavelength of the "extra-wide" lane carrier), using a combination of 
"pseudoranges" (or codes) Pew that share the same ionospheric dependence value and sign 
as the phase combinations. This is possible given the long wavelength resulting from the 
"extra-wide lane" combination. The values of. Lew and Pew can be determined by the 
following equations: 

/a Ji 

P.. = ^'^;^j^' =p*+a^J+M^ + E^^ + ... (2); 

equations in which Lx is the carrier phase observation (in units of length) with the 
frequency fx and with the wavelength Xx. The values of X 2ire shown in the attached 
TABLE I at the end of the present specification, it being understood that the numbers 1 
through 3 are associated with the three frequencies in the example described. TABLE I 
groups together a certain amount of associated data: frequencies, wavelengths, 
ionospheric dependence and maximum multipath errors, and various observations that are 
useful for a good understanding of the method according to the invention. This data is 
associated with the three frequency bands Z/ through Lj, respectively, and with the 
"extra-wide" and "wide" lane carrier phases {Lew and Zw), respectively. Px is the 
corresponding "pseudorange" observation. The multipath and observation errors for 
carrier phases and "pseudoranges" are marked nix, Mv, f^jc, and Z*^, respectively (the 
maximum multipath values and the typical measurement errors are also listed in Table I). 
The undifferentiated carrier phase ambiguity Xxbx, which contains instrumental delays 
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and, after double differentiation, becomes a multiplication integer value of the 
wavelength A^VANy, The ionospheric delay axl is proportional to the "Slant Total 
Electron Content," or "STEC," the integration of the ft^ee electron density along the beam 
radiated by a satellite normally being measured in TECU, as mentioned above. In 
equations (1) and (2), /?* represents a term that does not depend on the frequency 
(distance, clock errors, tropospheric refraction ,etc.). 

[0059] From equation (2), in which additional minor terms such as carrier phase 
wind'Up are not explicitly represented, it is possible to estimate the double difference of 
"extra-wide" lane ambiguity, hereinafter called VA/V^^^ , in a single observation "epoch" 
by subtracting the corresponding code, as expressed in the equation below: 

VAN^^ =1^ VA(4,, -P,J = VAA^,^ --?-(VAM,^ -h VA^,J + ... (3); 

[0060] When there are no substantial multipaths affecting the receivers (usually 
less than 3.7 meters), the error limit of such an evaluation is less than 0.5 cycles (as 
shown by the values provided in TABLE I), which makes it possible to instantaneously 
determine the ambiguity to an exact integer value. 

Step 2 : 

[0061] The combination "wide lane" ambiguity is estimated from the 
unambiguous "extra-wide lane" carrier phase obtained during the first step. Most of the 
time, the difference between the two is constituted by the "extra-wide lane" ambiguity 
and the differential ionospheric refraction (approximately 0.06 cycles/TECU with the 
working frequencies of the present example). The nondispersive terms are eliminated. 
The main problems here are the measurement error and the multipaths of the carrier phase 
signals. Although the typical values of the differential ionospheric refraction at middle 
latitudes and with baselines of less than 100 kilometers are just a few TECU, the use of 
the ionospheric correction can significantly increase the success rate at greater distances, 
and in more difficult ionospheric scenarios. 

[0062] The main details of this step will now be explained. 

[0063] Once the first "long wavelength" ambiguities VAA^^^^ are resolved using 

the corresponding "pseudoranges," a second ambiguity such as the "wide lane" 
combination Lw expressed by the following equation: 

A.. = A., -A)= ^'Y/f' =P* + KK + ccJ + m„^e„+... (4) ; 

Jl J 2 
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(an equation in which cpi and 92 are the phases with the frequencies Li and L2), can be 
estimated from VAN^^^ and from the corresponding difference of the carrier phases (see 
the corresponding constant values in TABLE I), in accordance with the equation: 

VAA^„, =^ VA(^,,, + m,^ -m^J +-i-(a^ J VA/ (5); 

In essence, in the presence of moderate multipaths, the error term corresponding to 
equation (5), which also includes the observation error, is normally less than 0.3 meters, 
i.e., less than 0.4 cycles. The error term remaining from the ambiguity resolution 
corresponds to ionospheric refraction at a value of 0.0580 cycles/TECU for the 
frequencies in the example (see TABLE I). This term may lower the success rate, but it is 
not a critical term at middle latitudes and at distances of less than a few hundred 
kilometers, for which the double differences of the "STEC" values, hereinafter called 
VA/, are generally less than 10 TECU. A detailed example will be given below in 
reference to the diagram in Fig. 5. 

[0064] It is in such a context that the use of a real-time ionospheric model 
improves the conditions for the success of "extra-wide lane" ambiguity resolution at 
lower latitudes and for scenarios with longer lines. 

Step 3 : 

[0065] The L/ phase ambiguity is derived from the difference between L/ and the 
unambiguous "extra-wide lane" obtained previously. In this step, the main problem is the 
corresponding differential ionospheric refiraction (approximately 1.9 cycles/TECU), 
which can produce errors of several cycles at middle latitudes. 

[0066] In the third step, an approach similar to the second step is implemented, 
but using the carrier phase differences between short and medium wavelengths instead of 
medium and long wavelengths, as is clear from the following equation (6): 

VAA^. =Y^ML, -A,+A,,7VJ = VAA^. _ J_ VA(^,, + -m,) + i-(a, -aJVA/ + ... ; 

an equation in which a/ and a,,> (generally ax) are ionospheric coefficients defined in 
TABLE I, in the example described. 

[0067] During the third step, the combination of the carrier phase measurement 
error and a moderate multipath introduces another error, typically less than 0.2 cycles (see 
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TABLE I). However, the critical problem here is the ionospheric refraction, which can 
also introduce errors of more than 0.5 cycles (-1.945 cycles/TECU) for short baselines. 

[0068] This fundamental limitation can be overcome by calculating real-time 
ionospheric corrections better than 0.26 TECU (i.e., 0.5 cycles x TECU/1.9475) in order 
to assure the correct integer ambiguity evaluation. These corrections are introduced in 
step 4, which will now be described. 

Step 4 : 

[0069] In order to overcome the problem that arises in the third step (i.e., in order 
to determine the ambiguities of the shortest wavelengths), a real-time model is 
determined in order to estimate the differential ionospheric refraction. This model is 
calculated from dual-frequency carrier phase data at fixed locations in a network of 
reference stations. It is assumed that a description of the ionospheric region passed 
through by the satellite-transmitted radio wave radiation is thus produced. The data of 
this model are combined in a known way with data derived from a simultaneous geodetic 
calculation, the latter advantageously being able to be performed by one of the fixed 
ground stations in the network, called a "master" station. The main advantage of this 
technique is that it makes it possible to estimate the differential refraction with a typical 
error of less than 0.25 TECU, even at distances of hundreds of kilometers from the 
nearest reference site and under varied ionospheric conditions. This precision typically 
involves errors of less than 0.5 cycles in the instantaneous determination of Lj ambiguity 
for medium and long distances. 

[0070] The ionospheric distribution of free electrons can be approximately 
determined by a grid of resolution volume units, or "voxels," wherein the electron density 
in an "Earth Centered Inertial" (or "ECI") system is presumed to be constant at a given 
moment. 

[0071] A typical setup of this type is described in Fig. lA. This Fig. lA 
schematically illustrates a meridian cross-section of voxels Voxyk^ (U j> k being coordinate 
subscripts for longitude, latitude and altitude, respectively), in which the ionospheric 
electron density distribution is broken down according to equation (7), so as to determine 
the data of a real-time ionospheric model. 

[0072] In this Fig. lA, the earth GT is shown in partial cross-section, surrounded 
by the ionospheric layer C/o/v, which has been arbitrary subdivided into two intermediate 
layers C,/ and C/2, respectively. The low altitude of the layer Cn is 60 km in the example 
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described, and its high altitude is 740 km. The low altitude of the layer Ci2 is 740 km, 
and its high altitude is 1420 km. The angular surfaces of the voxels Voxyk are 5x2 degrees 
each. 

[0073] By way of illustration, we have represented satellites SATj through SATn 
(three of which are shown in the figures), in which "GPS" transmitters, GPSei through 
GPSeh, respectively, are installed. These satellites SAT/ through SATn are in view of all or 
some of the ground stations, for example the three fixed reference stations ST/ through 
STm represented in Fig. I A, which include ''GPS" receivers GPSri through GPSro and 
integrated calculating means (not represented). 

[0074] Fig. IB schematically represents the architecture of a complete navigation 
system according to the invention, seen from the rover SUR. 

[0075] This system requires a computer, interfaces, and physical resources similar 
to those in certain currently deployed systems, such as EGNOS or WAAS. 

[0076] This system for implementing the method according to the invention 
includes a "GPS" receiver, SURgps, which receives the signals transmitted by satellites, 
for example the satellites SAT/ through SATn in Fig. lA. From these three- frequency 
band signals, the three steps of the "TCAR" method are performed as indicated above 
using integrated calculating means (not represented). 

[0077] It also receives the data signals firom the real-time ionospheric model 
determined by the reference station REF (the nearest of the stations in the network of 
fixed ground reference stations ST/ through STm of Fig. lA) from signals received from 
the satellites SAT/-GPSei through SATn-GPSEn using a "GPS" receiver, REFcps- The 
station REF transmits the data from the model to the rover SUR using a transmitter 
REFe/^ 

[0078] Geodetic configuration data, which can be determined in a single, so- 
called master station REF/^, are transmitted by its transmitter REF/^e and distributed to the 
rover SUR to be combined with the data fi*om the ionospheric model. As a rule, this 
station also has a "GPS" receiver REFMcps^ as do the other fixed ground reference 
stations. The geodetic configuration data are calculated simultaneously with the data of 
the geodetic model. 

[0079] In a way that is well known, the ionospheric determination can be obtained 
by resolving in real time, using a Kalman filter, the average electron density Ne of each 
illuminated cell and k {ij, and k being coordinate subscripts for three dimensions, as 
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defined above), treated like a random propagation process, and with a typical process 
noise of 10^ to lO'^ electrons/m3/ V/ic>wr . To give a nonlimiting example, a process using 
a Kalman filter is described in the book by G. J. BDERMAN entitled "Factorization 
Methods for Discrete Sequential Estimation," published in Mathematics in Science and 
Engineering, Vol. 128, Academic Press, New York, 1977. 

[0080] According to an important characteristic of the method of the invention, 
only the carrier phase data are used. Thus, "pseudorange" code noise and multipaths are 
avoided. The carrier phase polarizations 5/ (constant in each continuous arc of carrier 
phase data for each satellite-receiver pair) are estimated simultaneously, as random 
variables (which become random white noise processes when a cycle slip occurs). In the 
filter, the polarizations are decorrelated in real time fi-om the electron density values, 
inasmuch as the satellite geometry changes and the variations of the two types of 
unknowns have become smaller, as shown in equation (7) below, which represents the 
model for given ionospheric information between a "GNSS" satellite-receiver pair, Li and 
L2 being the carrier phases expressed in units of length, Li = Zy - L2, and Ne being the 
electron density). Li is given by the following equation: 

L, = STEC +B, = N^cll + = S Z S (^J.-../.* ^^/.y.* + ^/ 

an equation wherein the values "REC" and "SAT" relate to a GPS receiver and a satellite, 
respectively, of a "receiver-satellite" pair, dl is an elementary distance on the line of sight 
linking this pair, and As^. ^ j^. is an elementary coordinate surface of the STEC, and k 

being coordinate subscripts for three dimensions, as defined above. 

[0081] This approach is particularly suitable for detecting the local characteristics 
of the electron density distribution, and the use of two layers with ground data fi-om a 
Global Navigation Satellite System (abbreviated "GNSS") instead of a single layer as in 
the usual methods, significantly reduces the risk of a bad modeling of the electron 
content. 

[0082] In the case of "WADGNSS" ("Wide Area Differential GNSS") networks, 
it is possible, from the content of these real-time corrections of the slant total electron 
content, or "STEC" obtained through equation (7): 

[0083] - first, to form the station-satellite double difference VASTEC , with an 
error of less than 1 TECU, and to obtain a second (i.e., "wide lane") ambiguity in the 
reference stations; and 
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[0084] - secondly, to interpolate, in the rover's receiver, the unambiguous value 
Ly, i.e., a very precise value of VASTEC at a level of a few hundredths of a TECU, a 
value obtained after the resolution of the ambiguities. 

[0085] If the interpolated value is better than 0.26 TECU, then the rover's receiver 
can resolve both ambiguities in real time. 

[0086] This method is intrinsically similar to the so-called "WARTK" method, 
published in the aforementioned article by Hemandez-Pajares et al. (2000). The results 
obtained so far with the "WARTK" technique, in various experiments, are summarized in 
another article by Hemandez-Pajares et al., entitled "Tomographic Modeling of GNSS 
Ionospheric Corrections: Assessment and Real-Time Applications," published in "ION 
GPS 2001," September, 2001. It would be worthwhile to refer to these two articles for a 
more detailed description of this method. 

[0087] In the kinematic utilization of this method, one of the most important 
limitations is the existence of local ionospheric irregularities, such as propagating 
ionospheric disturbances, or "TID" (for "Traveling Ionospheric Disturbances"), which 
can produce incorrect results when using a linear interpolation of the ionospheric 
corrections between the reference stations in the network. Performance can be improved 
by incorporating dual-frequency data from the rover. 

[0088] Furthermore, in the case of the calculation in the reference station, long 
distances and strong electron content gradients can also limit the performance of the 
technique. For this reason, according to a characteristic of the method of the invention, 
an additional procedure was developed in order to include the case of permanent and/or 
reference stations separated by thousands of kilometers, using a so-called smoothed "wide 
lane" code, so as to improve the ionospheric model and consequently create an aid for 
real-time ambiguity determination. Such a procedure is described, for example, in the 
article by Hemandez-Pajares et al. entitled "Improving the Real-Time Ionospheric 
Determination from GPS Sites at Very Long Distances Over the Equator," published in 
Journal of Geophysical Research, 2002. 

[0089] However, the methods used in the prior art to resolve, in real time, the 
ambiguity in the receiver of a rover using a dual-frequency "GPS" system requires a large 
number of geodetic calculations in order to guarantee a relatively good position 
determination, on the level of twenty or so centimeters, to resolve free-floating 
ionospheric polarization combinations, and hence to be able (using the aforementioned 
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"WARTK" method) to resolve the ambiguities in real time, after a convergence time on 
the order of a few minutes. 

[0090] The invention makes it possible to overcome the deficiencies of the prior 
art, particularly as revealed by the "WARTK" method. It makes it possible to improve 
the instantaneous methods for resolving three-carrier phase ambiguities at medium and 
long distances between reference stations (ranging from tens to hundreds of kilometers), 
and with a minimal number of geodetic calculations. The main improvement lies in the 
third step of the aforementioned "TCAR" method, using real-time ionospheric corrections 
supplied by an ionospheric model calculated continuously in the reference stations. 

[0091] An additional improvement of the aforementioned step is also 
implemented in the method of the invention. This improvement consists in the use of 
three different pseudorange codes, hereinafter called Py, P^, Pj, instead of a single code 
element, which makes it possible to reduce the multipath effect in the resolution of 
"extra-wide lane" ambiguity. 

[0092] In addition, in a preferred embodiment, an integrity test is performed using 
the pseudorange, "wide lane," and Z/ frequency band codes to detect jumps in the 
ambiguity evaluation associated with an error in the ambiguity resolution of the longest 
wavelength. In essence, during the aforementioned second and third steps, it is easy to 
see from equations (5) and (6) that one cycle in the "extra wide lane" ambiguity error 
produces approximately eight cycles of "wide lane" error, and one cycle of "wide lane" 
error translates into approximately four cycles of Lj error , Li being the carrier with the 
shortest wavelength in the example described. Often, these jumps in the ambiguity 
evaluation are large enough to be detected and filtered using the corresponding codes. 

[0093] This approach according to the invention makes it possible to directly 
overcome the chief limitations that make a navigation with a typical error of less than a 
few centimeters impossible over long distances (more than 100 kilometers). 

[0094] TABLE II, located at the end of the present specification, summarizes the 
compared characteristics of the main methods for real-time ambiguity resolution in the 
prior art ("TCAR," "ITCAR," "WARTK") and according to the invention. 

[0095] To better illustrate the chief characteristics of the method according to the 
invention, we will now describe an exemplary experiment performed on the basis of 
several data sets supplied by modified signal generators installed aboard satellites in an 
existing system of the aforementioned "GNSS" type, with simulation of aircraft, roving 



16 



GRYN 221 (10505441) 

"surface" users (land vehicles, for example) or fixed sites. Various dynamic and 
ionospheric conditions were taken into consideration. 

[0096] By way of illustration and in order to provide a meaningful example, the 
real-time resolution of the ambiguity of a surface vehicle's receiver, hereinafter called a 
rover SURy was studied in detail relative to a fixed reference station, hereinafter called 
REF, separated fi:-om the rover SUR by approximately 129 km, in various scenarios. In 
the example described, two "GPS" carriers (1575.42 and 1227.60 MHz, respectively), and 
the 24-channel carrier of the "GLONASS" system, with a frequency of 1615.50 MHz, 
constitute the three frequencies adopted by a simulator of the "GNSS" system, for four 
satellites in view for 20 minutes at IHz, using the twelve available channels of a validated 
receiver of the type known as "AGGA" (for "Advanced GPS/GLONASS ASIC"). 
"AGGA" is a digital integrated circuit that allows high-speed digital signal processing for 
space applications such as radio occultations of "DNSS" system signals and low-earth 
orbit determinations. This circuit is based on an "ASIC" (or "Application Specific 
Integrated Circuit") standard, (an ATMEL T7905E component in the example described). 
"GLONASS" ("Global Orbiting Navigation Satellite System") is a Russian satellite 
navigation system that is comparable to the "GPS" system. 

[0097] Two basic data sets were considered in the example described, hereinafter 
referenced "P5-M0," a set associated with maximum signal power and no multipaths, 
which constitutes an ideal case, and "P3-M1," a set associated with average power and 
multipaths. 

[0098] In addition to the reference station receivers simulated in the data sets, 
three more stations existing in the "IGS"(for "International GPS Service") network, at 
distances of more than 200 km, were added just for the ionospheric calculation, which 
made it possible to emulate a more realistic situation comprising a larger network of fixed 
sites. 

[0099] We will now show the results of a detailed study for the most difficult 
case, i.e., the longest baseline of the "rover SUR - station REF"' pair, or approximately 
129 km. The reference station REF is one of the stations ST/ through STm in Fig. 1 A. 

[00100] Fig. 2 illustrates this example and represents the baseline and the 

positions of the station REF and the rover SUR relative to coordinate axes, ellipsoidal 
latitude and longitude, graduated in degrees. 
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[00101] Fig. 3 is an enlargement (zoom) of Fig. 2 showing the trajectory of 

the rover SUR, again relative to coordinate axes, ellipsoidal latitude and longitude, 
graduated in degrees. 

[00102] The first aspect of the results is referenced in terms of the real-time 

performance of the ionospheric filter. Examples of estimated delays from the real-time 
ionospheric model, referenced Tm in the graph, as compared to the real ionospheric 
delays, referenced Vv, are represented in Fig. 4. These values correspond to 
measurements from an experiment performed to show the capabilities of the method 
according to the invention. 

[00103] Also, parameters that are more important for precise navigation, 

the values and the differences (i.e., the errors) between the estimated double differences 

of "STEC", VA/ for the rover SUR, and their real values, VA/ , are shown in Figs. 5 and 
6 (on the ordinate, expressed in TECU; the abscissa representing delays in seconds). 

[00104] More precisely, Fig. 5 illustrates the double difference of "STEC," 

marked ''ddSTEC,'" estimated in real time with an ionospheric model according to one of 
the important aspects of the invention, compared to the true values for the receiver of the 
rover SUR, relative to the farthest receiver, or REF, (at a distance of approximately 129 
km). 

[00105] Fig, 6 illustrates the error in the determination by the real-time 

ionospheric model of the double difference of "STEC," again marked ''ddSTECr The real 
values of the ionospheric model are compared to two ionospheric thresholds, 
''Threshokir (+ 0.26 TECU) and "Threshold^ (- 0.26 TECU), making it possible to 
determine the third ambiguity VAA'^, . 

[00106] More than 92% of the ionospheric determinations have errors 

below the threshold limit of ± 0.26 TECU. They are therefore precise enough to make it 
possible to resolve the three ambiguities in the absence of multipaths, and to not take the 
measurement error into account. Most of the remaining 8% of evaluations having errors 
of more than 0.26 TECU came from a satellite that was observed at low altitude and in a 
southerly direction, the direction for which the ionospheric gradients are highest. These 
results were slightly improved by adding to the four reference stations a fifth station 
performing ionospheric correction calculations for the southerly direction. 

[00107] Once the real-time ionospheric corrections have been calculated 

and transmitted from the reference network, the steps of the method according to the 
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invention can be performed in the rover SUR in the manner described, in order to estimate 
and instantaneously resolve (in "single epoch" mode) the three ambiguities per beam. 

[00108] A summary of the main results of the ambiguity resolution, for the 

aforementioned data sets "P5-M0" and "PS-MI" respectively, is presented in the form of 
two TABLES, III and IV, located at the end of the present specification. The success rate 
of the three steps of the "TCAR" method are indicated in the following three cases: 

[00109] (a) without ionospheric corrections; 

[00110] (b) with corresponding ionospheric corrections fi*om a 

Klobuchar ionospheric correction model, transmitted by the current "GPS" system; 

[00111] (c) with corrections from a real-time ionospheric model, 

according to the method of the invention. 

[00112] More precisely, in the example described in connection with 

TABLE III, the success rate (in %) is relative to a total of 3834 on-the-fly resolution tests 
for the "extra-wide lane," "wide lane" and Li band fi-equency ambiguities (respectively 
VAA^^^^ , VAA'^^, and VAN^ ) for the receiver of the rover SUR, relative to the farthest fixed 

site REF (at a distance of approximately 129 kilometers), and with the aforementioned 
ideal data set P5-M0 (maximum reception power and no multipaths). 

[00113] TABLE IV is similar to TABLE III, but for the aforementioned 

data set P3-M1 (an average level of reception power and multipaths). The figures in 
parentheses indicate the success rate relative to the total number of observations. 

[00114] It is easy to see, when looking at TABLES III and IV, that the 

method of the invention, which combines the three steps of the "TCAR" method with a 
step for applying a precise real-time ionospheric model, makes it possible to considerably 
increase the success rate of instantaneous ("single-epoch" ) ambiguity determination, not 
only in an ideal scenario (data set P5-M0: firom 0 to 92%), but also in a more difficult 
scenario (data set "P3-M1") with multipaths (35% even using a "Klobuchar model" type 
of correction, whereas the success rate reaches 92% for VAN^ , with the corrections fi'om 
the ionospheric model produced according to the method of the invention). 

[00115] A look at TABLE IV shows that one of the important problems 

with the data set "PS-MI" affected by multipaths, is the lack of integrity in the ambiguity 
evaluations according to the "TCAR" method, VAA^^^^. and VAA^„,, with approximately 
10% and 4% incorrect instantaneous determinations, respectively. 
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[00116] To increase the integrity, i.e. to reduce the possibility of factoring 

in incorrect ambiguities, in this case in the presence of a carrier phase, and to code per 
multipath, it is possible to use the corresponding pseudorange, "wide lane" and Li codes 
to try to detect potential ambiguity determination errors in the ambiguity of the previous 
longest wavelength in the "TCAR" approach ("extra wide lane" and "wide lane" 
ambiguities, respectively), errors which are amplified by a factor of 9 and 4 wavelengths 
in the "extra wide lane" and "wide lane," respectively, as shown in equations (5) and (6). 

[001 17] A summary of the corresponding results is presented in the form of 

two TABLES, V and VI, located at the end of the present specification, respectively using 
the pseudoranges and the smoothed pseudoranges to filter such important errors. 

[00118] TABLE V is similar to TABLE II, but here the ambiguity 

resolution success rate is calculated after having been subjected to an integrity test, i.e., a 
comparison with the pseudorange-derived ambiguity, in order to filter the large jumps in 
"wide lane" and Zy ambiguity, these jumps being due to previous, larger wavelength 
ambiguity errors. The last column indicates the availability, i.e., the percentage of the 
3834 observations that pass the preceding integrity test with the pseudoranges. The figure 
in parentheses "*" indicates the percentages of ambiguities calculated after passage of the 
integrity test. As before, the figures in parenthesis indicate the success rate relative to the 
total number of observations. 

[00119] TABLE VI is similar to TABLE III, but smoothed pseudoranges 

are used in place of "pure" pseudoranges to perform the integrity test. 

[00120] It is easy to see that the method of the invention provides a 

significant improvement of approximately 20% in integrity (from 79% in TABLE IV to 
91%) in TABLE VI), with a relatively small reduction (16%)) in availability (which goes 
from 100%o to 84%), by using an integrity test with smoothed code. If instantaneous 
("single epoch") codes are used instead, the integrity is further improved (84%), and with 
greater availability (90%). When using a known "GPS" transmission model of the 
Klobuchar type in place of the corrections from the real-time ionospheric model 
according to the method of the invention, there is a substantial deterioration of the results, 
accompanied by a near-total lack of availability. 

[00121] The available data sets for the experiment just described, despite 

the fact that they were suitable for obtaining an instantaneous ambiguity resolution, were 
nevertheless very limited in determining an instantaneous ("single-epoch") position. In 
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fact, for the experiment, there were only four available satellites in view for 20 minutes of 
data, so the instantaneous navigation solution for the rover SUR had to be calculated 
using the data set "P5-M0" without multipaths, and with a very high available signal-to- 
noise ratio, in order to have the use of the minimum number of four satellites for 
navigation. Moreover, the tropospheric delay, the evaluation of which would have 
required at least a fifth satellite, was omitted. Furthermore, the first interval, from 
approximately 7200 to 7500 seconds, was skipped over in the determination of the 
positioning with respect to certain carrier phase cutoff configurations, potentially due to 
measurement problems in the receiver during this period. On the other hand, fi-om the 
point of view of the ionospheric corrections and the corresponding success rate of the 
ambiguity resolution, the results are practically the same as those obtained for the data set 
"P5-M0" (see TABLE III). 

[00122] These limitations are due to the use of data sets supplied by 

existing systems, which have been subject to only a minimum number of adaptations. 
With a more realistic data set, for example supplied by six or more satellites, these 
important limitations on the determination of positioning will disappear. It will in fact be 
possible to obtain a better geometry (a geometric "DOP" or "Dilution of Precision" 
factor) as well as the capability to detect and filter satellites (typically none, or sometimes 
one) having ambiguity errors. This will be able to be done using a navigation filter or 
autonomous receiver surveillance algorithms. 

[00123] The results corresponding to the four satellites contained in the data 

set P5-M0 are indicated in Figs. 7A through 7D, which show the residuals prior to 
adjustment (Fig. 7A), as well as the East (Fig. 7B), North (7C) and vertical (7D) 
components of the instantaneous positioning error after implementing the method 
according to the invention, for a rover SUR approximately 129 km away fi-om a reference 
station REF, The ordinates of the graphs are graduated in meters and the abscissas 
(delays) in seconds. 

[00124] For comparison. Figs. 8 A through 8D represent corresponding 

data, but for real ambiguities. 

[00125] The chief characteristic of the residuals prior to adjustment VAL^ 

(Fig. 7A) are errors of approximately 10 cm associated with an error of 1 cycle in L/ 
(approximately 8%: see TABLE III), which produce ("single-epoch") navigation errors 
similar to jumps. They are also represented in Figs. 7B through 7D. They are amplified 
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by the corresponding geometric "dilution of precision" factor ("DOP"). These incorrect 
evaluations of Li ambiguity typically affect one double difference out of each three that 
are available per epoch (the minimum number for determining a positioning), so this error 
affects the positioning three times (approximately 24% of the epochs, with 76% of the 
errors in the three dimensions under 5 cm and 100% under 21 centimeters). Although the 
distribution of the errors is not Gaussian, the resulting effective three-dimensional values 
("RMS") are 7, 3.5 and 2 centimeters, for the East, North and vertical components, 
respectively. 

[00126] To summarize these results, "RMS" errors of 1, 1 and 2 cm are 

obtained for the East, North and vertical components, respectively, v^hen the ambiguity is 
correctly resolved (approximately 92% of the tests and 77% of the epochs), and "RMS" 
errors of 3, 5 and 2 cm occur for the East, North and vertical components, respectively, 
when including the epochs with an incorrect ambiguity resolution. 

[00127] Furthermore, the observations reduced by residual calculations, 

prior to adjustment, of the carrier phase double difference and the instantaneous ("single- 
epoch") navigation errors are represented in Figs. 8A through 8D. A significant 
amplification of the carrier phase noise may be seen, particularly in the vertical 
component, producing an error of about 5 cm. This trajectory (with a three-dimensional 
"RMS" value of 3 cm, 1 .1 cm and 2 cm for the North, East and vertical components, and 
levels of 95% at approximately 2, 2 and 4 cm, respectively) very clearly represents the 
solution that may be obtained using a real-time filter instead of the instantaneous ("single 
epoch") mode solution, after correct resolution of the ambiguities during the first epochs. 

[00128] In order to characterize the impact of latency times in the 

ionospheric correction (due for example to potential problems in the communications), 
delays of 1 to 30 seconds (30 epochs) were considered for the ionospheric correction 
calculations in the network of fixed stations. For each of these delays, the success rate for 
a "ddSTEC" (in TECU) was calculated so as to obtain an absolute precision of 0.26 
TECU for a rover at a distance of 129 km. This rate is shown in Fig. 9 as a fimction of 
the latency time (in seconds). It is easy to see that the effect of the latency time is 
negligible up to 30 s, with a success rate of more than 90% achieved for all the latency 
times. 

[00129] In Fig. 10, longer latency times are considered, in this case for the 

reference network and rover corrections simultaneously, so that the ionospheric pierce 
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points remain close together. It is clear that after 5 minutes, the success rate decreases 
from 90% to 85%. After 10 minutes, the success rate drops fiirther, reaching 75%. These 
values can fall more sharply in scenarios wherein the variation in the differential 
ionospheric delay is even higher. 

[00130] It is clear that at a middle latitude, latency time does not pose an 

important problem for the method according to the invention, which can support typical 
latency times of up to about 5 minutes. 

[00131] Through the reading of the above, it is easy to see that the 

invention achieves the objects set forth. 

[00132] The method according to the invention makes it possible, by 

applying real-time ionospheric corrections using a tomographic model of the ionosphere, 
accompanied by a minimal number of geodetic calculations, to obtain a complete 
instantaneous ambiguity resolution. It also allows navigation with a precision on the 
order of one centimeter, at distances of more than 100 kilometers from reference stations. 

[00133] The preferred application of the invention will be in fixture satellite 

navigation systems such as "GALLILEO" or "Modernized GPS," which implement the 
transmission of carriers in three different frequency bands. 

[00134] According to a chief characteristic of the method according to the 

invention, this capability for instantaneous navigation, with a precision better than one 
decimeter, is due to the fact that a three-carrier phase ambiguity resolution method of the 
"TCAR" type, which is not very complex, is combined with a determination of a real- 
time ionospheric model in order to make it possible to apply ionospheric corrections. 

[00135] The exemplary experiments described in the present specification 

show a significant improvement of the results obtained by the methods of the prior art, 
under equivalent conditions: 

[00136] - The instantaneous ("single-epoch") success rate of approximately 

60% or less is improved, increasing to approximately 90% at distances of more than 100 
km from the nearest fixed ground reference station, and under conditions in which it is 
difficult to produce an ionospheric model (at noon at the solar maximum); 

[00137] - This high efficiency is maintained with latency times of 

approximately 5 minutes in the ionospheric corrections; 

[00138] - The corresponding navigation, despite the limited number of four 

satellites available, for the experimental conditions maintained by way of example, in the 
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data sets, can be achieved instantaneously, in "single-epoch" mode, and makes it possible 

to obtain three-dimensional "RMS" values of 3, 5 and 2 centimeters for the East, North 

and vertical components, respectively; 

[00139] - In a preferred embodiment, a real-time ambiguity integrity test 

made it possible, insofar as a simple pseudorange multipath reduction approach was 

incorporated into the method, to improve performance; and 

[00140] - This good performance was also proven in difficult scenarios: 

[00141] - Low latitude (35 degrees instead of 48 degrees), with higher 

ionospheric values. 

[00142] - Extreme conditions with respect to tropical stations beneath the 

"Northern Equatorial Anomaly" of the ionosphere, where the highest gradients are 
typically produced. Li the latter case, it is necessary for the user to produce its own 
ionospheric filter feed and combine it with the corrections from the reference networks 
and its own ionospheric observations. 

[00143] - High dynamics, making it possible to provide precise ionospheric 

corrections usable by an aircraft in flight at approximately 140 km from the nearest 
reference station. 

[00144] It should be clear, however, that the invention is not limited to just 

the exemplary embodiments explicitly described, particularly in connection with Figs. 1 
through 8. 

[00145] Finally, the numerical examples have been provided merely as an 

illustration, and do not constitute a limitation of any kind on the scope of the invention. 
They stem from a simple technological choice within the capability of one skilled in the 
art. 
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TAjBLEI 





X=l 


X=2 


X=3 


X=ew 


X=w 


Frequency 

(fx, MHz) 


1575.42 


1227.6 


1615.5 


40.08 


347.82 


Wavelength 

ttv, m) 


0.1903 


0.2442 


0.1856 


7.4799 


0.8619 


Phase measurement 
error 

fe, m) 


=0.002 


=0.002 


=0.002 


=0.1 


=0.01 


Code measurement 
error 

{Ex, m) 


=3 


=3 


=3 


=2.1 


=2.1 


Maximum phase 
multipath 

(«Jx, m) 


«0.05 
(<0.01) 


«0.06 
(<0.01) 


«0.05 
(<0.01) 


«2 
(<0.5) 


«0.2 
(<0.05) 


Maximum code 
multipath 

(Mv, m) 


«450 
(<10) 


«450 
(<10) 


«450 
(<10) 


«450 
(<10) 


«450 
(<10) 


Ionospheric 
coefficient 

(a^r, m/TECU) 


-0.1623 


-0.2673 


-0.1543 


-0.2083 


-0.2031 



TABLE II 





Advantages 


Disadvantages 


"TCAR" 


Low computational load. 


Ionospheric error, seriously 
limiting the resolution of the 

third ambiguity 


"ITCAR" 


Results improved by integrating 
TCAR into a navigation filter. 


The ionospheric delay continues 
to limit the resolution of the 
third ambiguity. 


"WARTK" 


A precise real-time modeling; 
allows precise navigation at one 
hundred kilometers from the 
nearest reference station. 


High computational load: a need 
to calculate a first ionosphere- 
free solution for the rover and a 
significant convergence time are 
necessary. 


Method 
according 

to the 
invention 


Low computational load and 
precise real-time modeling of the 
ionosphere, allowing precise 
"single-epoch" navigation at 
distances of up to 100 km from a 
reference station. 
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TABLE III 



P5-M0/SUR2-REF5 
(;=il29 km) 


Success rate 


Success rate 

VAN 


Success rate 

VAN, 


Without ionospheric 
corrections 


100% 


100% 


0% 


Klobuchar corrections 


100% 


100% 


33% 


Real-time ionospheric 
corrections 


100% 


100% 


92% 


TABLE IV 


P5-M0/SUR2-REF5 
(^^129 km) 


Success rate 

VAN 

ew 


Success rate 

VAN^ 


Success rate 

VAN, 


Without ionospheric 
corrections 


90% 


95% (86)% 


3% (2) 


Klobuchar corrections 


90% 


95% (87)% 


35% (31) 


Real-time ionospheric 
corrections 


90% 


95% (86)% 


92% (79) 



TABLE V 



P3-M1/SUR2-REF5 
(^129 km) 


Success 
rate 

VAN 

-* 


Success 
rate 


Success 
rate 

VAA^,* 


Availability 


Without ionospheric 
corrections 


99% 


95% (94) 


0(2)% 


38% 


Klobuchar corrections 


99% 


97% (96) 


33% (32) 


38% 


Real-time ionospheric 
corrections 


99% 


96% (92) 


91% (84) 


90% 



TABLE VI 



P3-M1/SUR2-REF5 
{^129 km) 


Success 
rate 

VAA^ * 


Success 
rate 


Success 
rate 

VAN, 


Availability 


Without ionospheric 
corrections 


100% 


100% (100) 


0% (0) 


0.002% 


Klobuchar corrections 


100% 


100% (100) 


37% (37) 


0.002% 


Real-time ionospheric 
corrections 


100% 


100% (100) 


91% (91) 


84% 



26 



